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In the circuit quantum electrodynamics architecture, both the resonance frequency and the cou¬ 
pling of superconducting qubits to microwave field modes can be controlled via external electric and 
magnetic fields to explore qubit - photon dynamics in a wide parameter range. Here, we experi¬ 
mentally demonstrate and analyze a scheme for tuning the coupling between a transmon qubit and 
a microwave resonator using a single coherent drive tone. We treat the transmon as a three-level 
system with the qubit subspace defined by the ground and the second excited states. If the drive fre¬ 
quency matches the difference between the resonator and the qubit frequency, a Jaynes-Cummings 
type interaction is induced, which is tunable both in amplitude and phase. We show that coupling 
strengths of about 10 MHz can be achieved in our setup, limited only by the anharmonicity of 
the transmon qubit. This scheme has been successfully used to generate microwave photons with 
controlled temporal shape [Pechal et ai, Phys. Rev. X 4, 041010 (2014)] and can be directly imple¬ 
mented with superconducting quantum devices featuring larger anharmonicity for higher coupling 
strengths. 

PACS numbers: 42.50.Ct, 03.67.Lx, 42.50.Pq, 85.35.Gv 


I. INTRODUCTION 


The strength of the interaction between an atom and 
the electromagnetic field is determined by the dipole mo¬ 
ment of the atom and the mode volume of the electro¬ 
magnetic field. Although in free space the interaction 
is typically weak, it can be enhanced by confining the 
field to a small volume in a cavity quantum electrody¬ 
namics (QED) setting. This field has seen tremendous 
progress m and has diversified from the traditional 
setting with real atoms to solid state realizations us¬ 
ing nanoscale electronic devices such as quantum dots 
Hig or superconducting circuits EHH] as artificial atoms. 
However, in most of the solid-state settings the coupling 
strength between the atom and the cavity modes is fixed 
by the geometry of the device and the position of the ar- 
tihcial atom in the cavity, neither of which can be mod¬ 
ified in situ. Although in recent years, superconducting 
circuit devices which allow in situ access to the ampli¬ 
tude of the qubit-resonator coupling have been realized 
[SHU], a scheme for controlling the phase of this cou¬ 
pling has only recently been demonstrated in [13) . Such 
a scheme is expected to be useful in a variety of settings, 
such as quantum gate operations [T3], creating shaped 
photons for quantum networks [131111 dl], measuring the 
vacuum state of a cavity exploring vacuum-induced 
Berry phases [T5], enabling the controlled coupling of a 
single or multiple qubits to multiple resonator modes m- 
121) . or engineering quantum reservoirs |22) . 
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The amplitude and phase tunability of the qubit- 
resonator coupling strength can be achieved by a two- 
photon process, known as a cavity-assisted Raman pro¬ 
cess, which has been extensively studied for A-systems 
P5H?7] and ladder-type (S) systems [H]- Cavity-assisted 
Raman processes employ an external coherent drive with 
a time-dependent amplitude fl cos {ujdt -I- 4>) to induce an 
effective coupling g between the qubit and the resonator 
degrees of freedom. The experimental access to the am¬ 
plitude (11) and the phase {(p) of the external drive en¬ 
ables in situ amplitude and phase tunability of this effec¬ 
tive qubit-resonator coupling. 

Cavity-assisted Raman processes can be readily 
applied to superconducting circuit elements of the 
transmon-type [13) which are in wide-spread use be¬ 
cause of their excellent coherence properties and the rel¬ 
ative simplicity of their fabrication. These circuit ele¬ 
ments realize an anharmonic oscillator system, i.e., a sys¬ 
tem in which transitions are allowed only between neigh¬ 
boring states and the transition frequencies differ from 
each other by multiples of a small negative parameter a 
which characterizes the anharmonicity. In experiments 
using the circuit QED architecture, not only the tran¬ 
sition between ground \g) and the hrst excited |e) state 
at frequency ujge, but also transitions between higher ly¬ 
ing energy levels can easily be addressed [50] and com¬ 
plex quantum states can be realized m- In particular, 
the second excited state |/), which is separated from |e) 
by uief = ojge + a, has been used widely for quantum 
gates [551135) . and plays an important role in our imple¬ 
mentation of the cavity-assisted Raman processes in a 
circuit QED setting. 

In our experiments, we investigate the tunability of a 
cavity-assisted Raman process induced coupling between 
a microwave resonator and a transmon device whose 
qubit states are dehned as the ground and second ex- 
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cited states. We demonstrate the amplitude tunability of 
the transmon-resonator coupling, and analyze the effects 
of the small anharmonicity of the transmon on the rate 
and fidelity of the population exchange (swap) between 
the transmon and the resonator, thereby complement¬ 
ing our experiments in which shaped microwave photons 
have been created and analyzed m- 

The outline of the paper is as follows. In Section Ull we 
present spectroscopic measurements of the tunable cou¬ 
pling strength g using a transmon-type superconducting 
qubit. In Section mil we derive an analytical expres¬ 
sion for g using first order perturbation theory in the 
drive amplitude Q,. In Section IIVI we explain an iter¬ 
ative method for calculating the drive-induced ac Stark 
shift of the qubit levels. In Section la we use numeri¬ 
cal simulations and second order perturbation theory to 
analyze the fidelity of the excitation exchange between 
the transmon and the resonator induced by the tunable 
coupling. 


II. SPECTROSCOPIC MEASUREMENT OF 

THE TUNABLE TRANSMON-RESONATOR 
COUPLING 

Our system consists of a transmon-type supercon¬ 
ducting qubit with maximal Josephson coupling energy 
^max/^ _ ^ GHz and charging energy Ec/h = 

0.343 GHz. We operate the transmon at a transition 
frequency of = 8.103 GHz which is higher than 

the fundamental mode frequency a;r/27r = 7.126 GHz of 
the resonator, resulting in a positive transmon-resonator 
detuning A = Wge — a;,. = 27r x 0.977 GHz. With a 
coupling strength g/27r = 65 MHz between the trans¬ 
mon g-e transition and the fundamental resonator mode, 
the system is far in the dispersive regime (A 3> g). The 
transmon has an anharmonicity al2'K = — 0.376 GHz, 
and the frequency of the transition between the first and 
second excited state is We/ = Wge + a = 2 tt x 
7.727 GHz. The resonator decay rate k/2'k is measured 
to be 6.6 MHz. 

A tunable effective Jaynes-Gummings-like coupling be¬ 
tween the transmon state |/) and the resonator can be 
activated by applying a coherent microwave tone at a fre¬ 
quency close to the energy difference between the dressed 
states |/ 0 )d and \gl)D- Here, the subscript D denotes 
the dressing of the combined eigenstates |/ 0 ) = |/) 0 | 0 ) 
and I 5 I) = I 5 ) 0 | 1 ) by the transmon-resonator coupling 
g. The effective coupling leads to a coherent excitation 
exchange between the transmon and the resonator. The 
fidelity of this exchange is maximum when the transmon 
drive frequency is equal to wj], which is defined by the 
resonance condition 

= 2ujge + a — Ulr E ^fOgli^)- (1) 

Thus, is the angular frequency difference between 
\fO)D and \gl)D modified by the difference A/ogi(f^) be¬ 


tween the ac Stark shifts for \fO)D and \gl)D, which de¬ 
pends on the amplitude of the coherent drive H. 
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FIG. 1: (a) The resonator transmission as a function of 

transmon and resonator drive frequencies. The white curve 
is a fit to the data at the fixed transmon drive frequency 
LOdl^Tv = 8.75 GHz. (b) The normalized transmission data 
for the white fit curve in (a). The amplitude of the effective 
coupling strength g/2n extracted from this trace is 3.1 MHz. 
(c) The measured effective coupling g (diamonds) as a func¬ 
tion of coherent drive amplitude H, the perturbation theory 
prediction (thick red line) and the numerical simulations (thin 
green line). 


The strength of the effective transmon-resonator cou¬ 
pling g is measured by weakly probing the frequency- 
dependent transmission through the resonator [ 6 ] . In the 
strong-coupling regime, the transmission peak of the res¬ 
onator splits into two distinct peaks of equal width when 
the resonance condition in Eq. 0 is satisfied, as ob¬ 
served in the measurement data shown in Fig. 0 a-b). 
The frequency separation between the two maxima in 
transmission then equals 2gf2TT, i.e., twice the coupling 
strength between the states |/0 )_d and \gl)D- 

To extract the effective coupling we first identify the 
transmon drive frequency wj] at which the transmission 
curve is split into two peaks of identical width. We then 
fit the resonator transmission curve to the response func¬ 
tion of the transmon-resonator system given by 


S{ojd) = Al 


4g2(a;0)2 


* I7I w — (D^ 

{i I 7 I w — |k| uj — (D^) 


( 2 ) 


derived from the master equation of the coupled 
transmon-resonator system in a truncated basis |36j . In 
Eq. 0, 7 is the qubit decay rate, and 

We have measured the coupling strength g for in¬ 
creasing amplitude of the coherent drive strength H 
[see Fig. 0c)]. When the drive is weak, the effective cou¬ 
pling strength g increases linearly with the amplitude of 
the drive strength, in good agreement with a first-order 
perturbation theory calculation outlined in section [ml 
However, for drive amplitudes larger than approximately 
0.2 GHz, higher-order effects start to contribute signifi¬ 
cantly and the dependence of 5 on H becomes non-linear, 
making the g smaller than predicted by the linear model. 
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The measured strength of the tunable coupling shows 
good agreement with the numerical simulation [Fig.j^c), 
thin green line] and the analytical result from a perturba¬ 
tion theory calculation to first order in [Fig. [^c), thick 
red line], which are discussed in the following. All system 
parameters which are relevant to the calculations were ex¬ 
tracted from separate experiments. The drive amplitude 
seen by the transmon qubit was calibrated by fitting the 
n dependence of A/ogi- The calibration routine for the 
coherent drive strength is further discussed in Section HVl 


III. CALCULATING THE COUPLING 
STRENGTH g 

In a reference frame rotating at the frequency the 
Hamiltonian for the transmon coupled to a resonator 
mode can be written as a sum of an N-level Jaynes- 
Cummings term and a coherent drive term, H = Hjc + 
Hd. Here [29], 


H]c = + Sqb^b -I- —b^b^bb (3) 

+g{ab^ + a^b), 

and 

Hd= ^(e'-^fe + Ac.) , (4) 

where Sr = uj^ — uJd and 5q = cUg^ — cud as the resonator- 
drive and transmon-drive detunings, respectively, w® and 
Wgg denote the bare quantities associated with uJr and cuge 
defined above (Fig. [^. The operator a (a'^) is the anni¬ 
hilation (creation) operator for the resonator mode and 
b(b^) its analogue for the transmon b = ]g)(el+-\/ 2 Je) (/]-!- 
•\/3|/)(fi| -I- ... when treated as an anharmonic oscillator 
at frequency (35). The Jaynes-Cummings type inter¬ 
action m couples the states which have the same total 
number of excitations, while the coherent drive field with 
time-dependent amplitude n(t), frequency Wd, and phase 
(j) couples neighboring pairs of transmon states. In the 
following, we omit the time dependence of H(t) for no- 
tational clarity. The above Hamiltonian is valid in the 
transmon limit where E]/Ec ^ 1 - 

We consider the system in the dispersive regime, that 
is, A = ojge — ojr ^ g- To calculate the tunable effective 
coupling strength, we rewrite the Hamiltonian E[ in the 
eigenbasis of FIjc, and treat the coherent drive term Eld 
perturbatively, expanding the solution in powers of the 
small parameter ft. Note that \i,j)D (* = g,s.^ f,h,... 
and j = 0 ,1, 2 ,...) is used for the eigenstates of iFjc, 
and \i,j) for the bare qubit-resonator states. The drive 
Hamiltonian Eld can be written in the \i,j)D basis as 

ij,kl 



FIG. 2: Schematic energy level diagram of the system (a) in 
the laboratory frame and (b) in the rotating frame of the drive 
frequency iJd that compensates for the energy difference be¬ 
tween I/O) and |( 7 l). The effective tunable coupling g between 
the bare states |/0) and |grl) can be understood as a second 
order cavity-assisted Raman process, indicated by the dashed 
arrows. 


from the terms coupling the resonant states \g, Z-l-I )d and 
\f,l)D while the terms describing off-resonant transitions 
in Eq. (§ can be neglected. With this rotating-wave-type 
approximation, the drive Hamiltonian becomes 

Hd « Y.^\9J + i)D{f,l\D 

I 


To show that the effective coupling between the trans¬ 
mon and the resonator is indeed of the Jaynes-Cummings 
type, we note that the dressed states ]/, 1)d and \g, Z + 1 )d 
are given, up to first order in g, by 


\f,l)D 


Iff, I + l)n 


I r 7\ ff-\/2(/ -I-1) , , 


gV^ 

' A + 2a 




Iff, Z -I-1) + 


ff-\/r+T 


A 


e, /), 


( 6 ) 

(7) 


where the coupling g{ab^ -I- a^b) is considered as a per¬ 
turbation to the uncoupled transmon-resonator system 
(i.e., (ff/A)^ ^ 1). Using this approximation for the 

dressed states, we calculate the matrix element 


~gi ^ (ff, l+MoHdlf, 1)d ^ g^^'^\/ ^A(A''+ a) ^ 

which represents the coupling strength between dressed 
states I/, and |ff, / -I- 1)d- As expected, gi is tunable 
both in phase and amplitude because of its dependence 
on the complex drive strength 

Next, defining the qubit raising and lowering operators 


When the resonance condition in Eq. 0 is satisfied the 
main contribution to the evolution of the system comes 


^ = '^\f,l)D{g,l\D, (9) 

I I 
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and the dressed photon annihilation operator 


— X] \jJ)D{jJ + Mdi (10) 




the drive Hamiltonian in Eq. can be written in the 
Jaynes-Cummings form 


Hd « gd^ai+g*d(7i'’ 


( 11 ) 


with 5 = go Ki a/(\/2A{A + a)). The absolute 

value of the coupling g in Eq. ( [TT| ) descriloes the splitting 
observed in Fig. [^c). 

It is instructive to compare the tunable transmon¬ 
resonator coupling g to the effective coupling g/^ = 
gfle^'^/A between the two degenerate states of a A sys¬ 
tem obtained from the adiabatic elimination technique 
[2?1 E5l [551 IMl- Most notably, the coupling strength 
for the transmon-resonator system is lower than that 
of the A system, g < gA- This result follows from the 
opposite signs of perturbative contributions from \e,l) 
and |e, Z -I- 1) in Eq. (^). Physically, this effect arises 
from the destructive interference of the two second order 
transition paths that couple degenerate levels \fO)D and 
l5l)i> (Fig-[2|- Indeed, the only reason that the coupling 
g does not vanish completely is the anharmonicity a of 
the transmon qubit, which results in a difference in the 
magnitude of perturbative contributions from |e, Z) and 
|e, Z -I- 1). On the other hand, in a A system, there is 
only one transition path coupling the degenerate levels, 
and consequently there are no interference effects. The 
correspondence between the A and the transmon systems 
is easy to see in the limit of large anharmonicity, when 
g[a. ^ oo) = gA- 

As shown in Fig. Be), the first order perturbation the¬ 
ory gives a satisfactory approximation to the tunable 
transmon-resonator coupling for small drive amplitudes 
which satisfy (H/2(A-|-a)) ^ 1. However, the ap¬ 

proximation to g breaks down when this inequality is 
no longer satisfied, and the first order approximations to 
\fO)D and \gl)D in Eq. (^) lose their validity. In par¬ 
ticular, we observe in the experiment that the first order 
approximation starts to break down as H/27r is increased 
above approximately 0.2 GHz for A/27r = 0.979 GHz, 
i.e., (rZ/2(A-|-a)))^ « 0.025. 


transmon qubit we fit the observed Stark shift to a per¬ 
turbative expression for A fQgi , where all the parameters 
other than H can be determined from separate measure¬ 
ments. In the following, we discuss the resolvent method 
used to obtain such a perturbative expression for A/ogi- 

The resolvent method allows for a systematic approx¬ 
imation of Afogi for increasing orders of interaction in 
both g and H. The Hamiltonian H = Hq -|- Hi consists 
of a bare part Hq and an interaction part Hj and its 
resolvent is defined as [iO] : 


G{z) = 



( 12 ) 


where z is a complex variable. In particular, the eigenen- 
ergies of H are given by the poles of G{z). Specifically 
for our transmon-resonator system, we split the driven 
Jaynes-Gummings Hamiltonian given by Eqs. (§ and 0 
into 


Ho = Sro) a + Sqb^h (13) 

and 

Hj = g ab^+ + h.c., (14) 


where g and H are real numbers. With this notation, we 
can rewrite the resolvent as 



(15) 


If the states \(j)i) and \<j)j) are two degenerate eigenvectors 
of the bare Hamiltonian satisfying 
the resolvent operator restricted to the two-dimensional 
space spanned by them can be written in the final form 


G{z) = ---—, 

z- Ho - E(z) 


(16) 


IV. AC STARK SHIFT AND DRIVE POWER 
CALIBRATION 

To determine the conversion factor between the ap¬ 
plied drive power and the drive amplitude H seen by the 


where the operators Hq and E(z) are to be understood 
as acting on the two-dimensional space and the matrix 
elements Atki{z) = {(j)k\^{z)\(l)i) for k,l G {Z,j} are given 
by 
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Sfci(z) = 


/ . . -^m 

m^ 2 j 

+ ^ ' {4>k\Hl\(j)m) _ p _ p/ ('/'m' + ’ ’ ’ , 

,/■ ■ ^ ^ ^ra 

m 'm ’ -/-■> 'I " 


(17) 


which is the weighted sum of all transition paths coupling 
states \(j)k) and through intermediate bare eigenstates 
{\4’m)}m=i excluding \(j)i) and |(/)j). Equation (16) is ob¬ 


tained by inserting an appropriate number of copies of 
the identity I = \4>m){4>m\ into each term of the sum 

in Eq. (15), and by noticing that the identity 




= ('('fcl 



holds if fc, Z G {i,j}- 

Eq. (16) shows that Hq + S(z) can be interpreted as 
an effective Hamiltonian H^^{z) describing the evolution 
of the system in the two-dimensional subspace spanned 
by \(j)i) and The resonance condition ([T| ), which 

corresponds to the avoided crossing in Eig. [l] can now 
be expressed as 

^ - Hgigi{z) = z- i7/o/o(^) = 0> 


which is to be satisfied when the bare energies Efg and 
Egi are separated by a detuning AfQgi{fl) + Ajq. Here 
the constant term Ajc represents the H-independent 
renormalization of the bare energies due to the Jaynes- 
Cummings coupling g only, such that A/ogi(0) = 0. 
Thus, we can calculate the detuning Afogi{fl) by iter¬ 
atively solving the coupled equations 


Z- Egl- Eglgl{z, H) = 0 

2 : — Egl + Afogl{^l) + Ajc — T,fofo{z, ft) = 0, 

(19) 


with the initial value z = 0. Eitting the resulting expres¬ 
sion for Afogi{ft) to the experimental data (see Fig. 
provides the conversion factor k between the applied 
drive power and the drive amplitude ft seen by the trans- 
mon. This conversion factor is also used to compute the 
qubit drive strength value in Fig. 0c). giving good agree¬ 
ment between our measurement and perturbation theory 
calculation. 


V. NUMERICAL AND ANALYTICAL RESULTS 
ON TIME-DEPENDENT VACUUM RABI 
OSCILLATIONS 

The effective transmon-resonator coupling leads to a 
coherent exchange of excitations (swap) between the 


states \f,l)D and \g,l + 1)d- Achieving a high fidelity 
for this swap operation is crucial for further applications 
of the microwave-induced transmon-resonator coupling in 
the context of quantum computing. In this section, we 
first derive an analytical expression for the fidelity of the 
swap operation when the external drive is turned on and 
off instantaneously. Then we use numerical simulations 
to account for the effects arising from pulse profiles that 
vary slowly with respect to A. We conclude that the 
swap operation can be realized with very high fidelity for 
realistic pulse profiles. 

When the drive is turned on and off instantaneously, 
one can derive an analytical expression for the Rabi os¬ 
cillations in the population in \gl)r), given the initial 
state \fO)D- The population in \gl)D at time t, denoted 
Pgioit), is given by the modulus square of the overlap 
between the time evolved initial state |/0(t))D and the 
target state |gl)_D 

PgiD{t) = \{gl\D\im)D\\ (20) 

where |/0(t ))d = U{t)\f0) d with the unitary time evo- 


qubit drive power P [dBm] 


-00 02 4 6 8 10 12 14 16 18 



qubit drive Strength, fi/ 27 r [GHz] 


FIG. 3: Measured values of A/ogi for increasing drive ampli¬ 
tudes ft. A fit to the expression in Eq. | |19[ ) obtained from 
the resolvent method (solid line) where the pertnrbative se¬ 
ries is truncated at the the sixth order in ft. The fit gives 
the conversion factor between the applied drive power and 
the applied drive amplitude seen by the transmon. The scale 
on the top indicates the drive power at the signal generator 
where each tick stands for the drive power at which the data 
point is obtained. 
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lution operator U{t) = 

To express Pgioit) in a convenient form, we expand 
the initial and target state in the eigenbasis of the full 
Hamiltonian 


|/0)p=^a„|$„) 

n 

|gi)D = 


( 21 ) 

( 22 ) 


where the eigenstates satisfy = eni'hn), and 

the coefficients a„ and /3„ are defined as ($n|/0)D and 
(^n|ffl)_D, respectively. In this new basis, the time- 
dependent population in \gl)D is 

PglDit) = (23) 

n 

+2 X! \^*ml^rnall3n\ COS [(e„ - e„)t -h 9nm] , 


where Onm is the phase of {a!^f5rnCt*nl3n)■ 

Equation (23) is an exact and useful relation between 


the solution to the full Hamiltonian (which can be ap¬ 
proximated using perturbation theory) and the rate and 
fidelity of the population exchange between |/0)d and 
\9^)d- In particular, the maximum fidelity P of the 
swap operation between |/0 )d and \gP)D can be ex¬ 
tracted from Eq. (23) as = 4|a+a_/3+/3_| and its 


rate as 2g = e+ — e_, where the subscripts ± stand for 
the two eigenstates of which have the largest overlaps 
with |/0)_D and \gl)D- We denote these eigenstates with 
|$±). When the resonance condition in Eq. ([^ is satis¬ 
fied and there are no other states whose energies are close 
to that of \fO)D or |gl)_D, i.e., (II/2/(A-|-a))^ < 1, 

these eigenstates are simply the two polariton states 
|$±) = 1/^/2 {\f0)o ± | 51 )d)- In this weak drive limit, 
Eq. ( [^ implies that the population exchange between 
|/0)d and \gl)D occurs with unit fidelity P = 1. 

For realistic amplitudes and detunings a reduction in 
P is caused by the population leakage out of the ini¬ 
tial and target states. To calculate the fl-dependence 
of the fidelity, we use the second order corrections to the 
eigenstates |<I>±) induced by the perturbation Hd (see ap¬ 
pendix!^. As a result, the fidelity of the swap operation 
between[/0)£) and \gl)D is 


P=l- 


VL/2 
A + a 


(^/2n/2V (^/in/2\ 


(24) 

up to second order in H/A. Notice that the reduction in 
fidelity is caused by population leakage out of the states 
\fO)D and \gl)D to neighboring states. We observe that 
the simulated fidelity shown in Fig. for a drive pulse 
with zero rise time is in good agreement with the an¬ 
alytical calculation in Eq. (24). The simulation of the 


transmon dynamics in the absence of decoherence was 
performed by solving the Schrodinger equation with the 
Hamiltonian given by Eqs. (§ and ([^. The underlying 


Hilbert space was truncated to four energy levels of the 
transmon and four Fock states of the resonator. 

The calculation above is based on the assumption that 
the drive pulse is turned on an off instantaneously and 
the obtained fidelity is therefore valid in the limit of an 
ideal square pulse. If, however, the pulse is ramped up 
and down gradually, the process becomes adiabatic with 
respect to the off-resonant transitions. Simulation results 
presented in Fig. show that this effect improves the 
fidelity significantly. 

The shape of the drive pulse 11 (t) for the simulations is 
chosen such that the effective coupling g{t) has the form 

9{t) = ffmax for At < t < r - At, 

= ffmaxSin^(7rt/2At) for t < At, 

= ffmaxSin^(7r(T — t)/2At) for t>T— At, 

where T denotes the length of the pulse. At the rise time 
and the amplitude ^max is chosen such that g{t) dt = 
tt/2, resulting in a 7r-flip between the states |/0 )d and 
The frequency of the pulse is varied in time to 
exactly cancel the variations in the amplitude-dependent 
ac Stark shift. 

To generate the correct drive pulse for the simulation, 
we needed to calculate the ac Stark shift and the effec¬ 
tive coupling g with good accuracy. For this reason, we 
decided to use a numerical procedure based on diagonal- 
ization of the Hamiltonian which is described in more de¬ 
tail in Appendix While this method is more accurate 
than the analytical expression obtained using perturba¬ 
tion theory, the latter is a more time efficient way of 
determining the behaviour of ac Stark shift as a function 
of other system parameters (i.e. A, a, and g). 

With pulses generated by this method, the simulations 
show that the population swap between |/0)d and \gi)D 
can be realized with essentially unit fidelity (> 0.99999) 
when decoherence is neglected. In current state-of-the- 
art experiments, the achievable fidelities will therefore 
be limited mainly by coherence times of the transmons 
or imperfections in the generated drive pulses. 


VI. CONCLUSION 

In conclusion, we have experimentally demonstrated 
a cavity-assisted Raman process to realize a tunable 
transmon-resonator coupling in a superconducting circuit 
QED architecture. This effective coupling is induced by 
coherently driving a transmon qubit. Its maximal value 
is only determined by the maximum qubit drive power 
that can be applied to the qubit. The measured data is 
in very good agreement with a perturbative calculation 
of the coupling strength g between |/0 )d and \gl)D- Our 
calculations show that the strength of the Raman tran¬ 
sition is reduced in comparison to the one implemented 
in a A level system [SUET] by the destructive interfer¬ 
ence between the two second order transition paths which 
couple |/0 )d and \gl)D- 
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FIG. 4: (a) Dependence of the loss of fidelity 1 — 7^ on 

the length of the tt pulse and its rise time. The transmon¬ 
resonator detuning is set to A/27r = 0.979 GHz, the anhar- 
monicity is a/27r = — 0.376 GHz, and gr/27r = 65 MHz. 
The solid lines show the results of the numerical simulations. 


The theoretical prediction of Eq. (241 is shown in dashed red. 
(b)-(c) Time evolution of the populations pfo, Pgi in |/0 )_d, 
|51)d, and the leakage to other levels given by 1 — pgi — Pfo 
for a 50 ns drive pulse with rise times of (b) 0.2 ns and (c) 
5.0 ns. 


We also determined the fidelity of the swap opera¬ 
tion which utilized the cavity-assisted Raman process, by 
both numerical and analytical means. A second order de¬ 
generate perturbation calculation shows that the fidelity 
of the population exchange strongly depends on the pop¬ 


ulation leakage from |/0 )_d and |gl)_D to the closest lying 


states [see Eq. (24)], a result which is also supported by 
numerical simulations (Fig. [^. We found that the fi¬ 
delity of the swap operation is expected to be very close 
to unity in the absence of decoherence effect and in real¬ 
istic systems will most likely be limited by the coherence 
time of the transmon. 
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Appendix A: Numerical calculation of the ac Stark 

shift 

In this appendix, we describe the evolution of the 
transmon under the approximation of slowly varying 
drive pulses. We then define the ac Stark shift opera¬ 
tionally as the shift of the drive frequency required to 
implement a perfect swap between the states \f0)D and 
Finally, we describe a method for calculating the 
ac Stark shift numerically. 

If the drive amplitude 17 (t) is varied slowly compared 
to the energy separation of |/0)d and |gl)_D from other 
energy eigenstates, the system evolves adiabatically with 
respect to the off-resonant transitions. Therefore, if the 
system is initially prepared in a superposition of \f0)D 
and |(?1)_D, it remains at all times in the subspace 5(17) 
spanned by the two instantaneous eigenstates |<i)i(17)) 
and 14)2(17)) corresponding to |/0 )_d and \gl)D- 

To describe the evolution of the state vector |'I'(i)) in 
5(17), we first introduce a mapping which connects the 
subspaces 5(17) for different values of 17. Using this trans¬ 
formation, the action of the full Hamiltonian for any 17 is 
mapped to an effective Hamiltonian acting on 5 (0), that 
is, the subspace spanned by \f0)D and |gl)_D. Using this 
formalism, we calculate the Stark shift as the shift of the 
drive frequency needed to keep the effective Hamiltonian 
in the “resonant form” ga^ = ff(|/0)D(5l|D + H.c.). If 
the drive frequency is not adjusted, the effective Hamil¬ 
tonian contains a term proportional to which prevents 
us from realizing a perfect swap operation between the 
two states. 

We start by introducing a linear map Mq : 5(0) —>■ 
5(17) defined as 

Mq = lim P(17)P(17 - A17)... P(2A17)P(A17)P(0), 

AO—>0 

(Al) 

where P(17) = X]i=i 2 l‘^i(^))(‘^i(^)l is a projector onto 
5(17). This map represents a continuous series of pro¬ 
jections onto the subspaces S{x) for x varying from 0 to 
17. 

As the state vector |4'(7)) evolves adiabatically and 
therefore lies in 5(17(7)), we can write it in the form 

|T'(7)) = Mo(t)|V)(7)), 




















































where is some vector evolving in the subspace 5(0) 
spanned by |/0)d and \gl)D- In this way, we have re¬ 
duced the problem of finding the evolution of the state 
vector l'I'(t)) in a changing subspace 5(n(t)) to that of 
finding the evolution of the directly related vector \tp{t)) 
in a fixed subspace 5(0). 

Since |d>(t-|-dt)) lies in 5(f2(t-|-dt)), it can be expressed 
as 


|'l>(t -I- dt)) =P{n{t + df))|'I'(t -I- dt)) 

=P{fl(t + dt)) exp(—ii?(t)dt)|'P(t)), 

resulting in the following evolution equation for |'0(t)): 


Mfi(^t+dt)\'4’it d-dt)) — (A2) 

P(n(t + dt)) exp(-iP(t)dt)Mo(t)|'!/'(t)) 

Under the reasonable assumption that the subspace 5(U) 
changes smoothly with U, it can be shown that Mq pre¬ 
serves vector norms in 5(U). It follows that which 
is an infinite product of projectors analogous to Mq but 
in the reverse order, is the inverse of Mq. Hence, after 
dropping the pro jector P(r2(t -|- dt)) from the right-hand 
side of Eq. (|A2|) as well as from the product form of 


Mn{t+dt) (cf. Eq. (Al)) on its left-hand side, we multiply 


the equation by to obtain 


+ dt)) = exp(-iP {t)dt)M^f^^^ 


which we transform into the differential form 

^|V-(t))=-iMt(,)P(t)M^(,)|V'(t)). 

The evolution of the vector \ip{t)) is therefore governed 
by an effective Hamiltonian 


that the equal superposition states |<pi, 2 ) = (|/0)d ± 
\gl)D)/V2 are eig enstates of Hee{t) and therefore, by 
virtue of Eq. (A3), that MQ^f^\(pi^ 2 ) are eigenstates of 
P(t) which we have previously denoted by |fl>i, 2 (U)). In 
other words, 


|$,(H)) =Ma|vJ.). 


(AS) 


This equation can be solved for Mq. However, since 
our goal is to determine the ac Stark shift, we need an 
equation for the Hamiltonian instead. To get it, we trans¬ 
form Eq. (AS) into a differential form. By substituting 
U —>■ H -I- dU, we find the following relation between 
|$i(U-|-dH)) and |$i(H)): 

|$i(H -l- dU)) = P(H -|- dH)|d>i(H)), 

which, after multiplication by ($j(U + dU)!, we write in 
the form 


($,(H)|A|$^.(f^))=0. 


The associated initial condition follows from substitut¬ 


ing H = 0 into Eq. (AS), giving |$i,2(0)) = (|/0 )d ± 


\gl) £i)/-\/ 2. Therefore, the two equal superposition states 
have to be eigenstates of the non-driven Hamiltonian. 
This is by definition also true for \fO)D and \gl)D- The 
only way the two distinct pairs of vectors can be eigen¬ 
states at the same time is if the subspace they are span¬ 
ning is degenerate. This can be achieved by choosing 
the correct freqnency of the rotating frame, giving us a 
condition for the drive frequency at H = 0. 

For i = j, the differential equation above can be sat¬ 
isfied simply by choosing the correct phase of the eigen¬ 
states |<i)i_ 2 (H)). After expressing the derivative of the 
eigenstate in terms of the derivative of the Hamiltonian, 
the remaining equations for i ^ j are equivalent to 


Hesit) = (A3) 

acting on 5(0), resulting in the evolution operator 
P(tf, ti) = T exp /j/ (—iPeff(i)) dt. If we assume that the 
drive amplitude at the initial time ti and the final time 
tf is zero, we have Mq^i.') = MQ(tt) = T’(O) and hence 
|'0(ti)) = |dt(ti)) and |'0(tf)) = l'k(tf)). Then we can 
directly write down the evolution from |'I'(ti)) to |4'(tf)): 

|vl/(tf)) = p(tf,ti)|vl-(ti)). 


This result allows to determine the ac Stark shift as the 
amplitude-dependent shift of the drive frequency needed 
to realize a perfect swap operation between \fO)D and 
For this, the effective Hamiltonian has to have 
the form 

Peff(t) =Poffset(i)l+5(0(l/0)D(5lb+H.c.), (A4) 

where the overall energy shift Poffset(i) leading only to 
an overall phase shift is omitted since it is physically ir¬ 
relevant. This equation is equivalent to the requirement 




dH{n) 

dn 


|$2(fl)) =0. 


Now we consider that the Hamiltonian depends on 
H not only directly but also through the amplitude- 
dependent drive frequency uid{^). The total derivative 
with respect to U can then be expressed using the chain- 
rule, leading to the equation 


(<i>i(U,a;d)|^^^^|^|$2(fI,a;,)) + 

^(^($i(H,a.,)|^^£p^|<i>2(fI,a;d))=0. (A6) 


We can solve this differential equation for Wd(H) nu¬ 
merically to find the amplitude-dependent drive fre¬ 
quency which yields an effective Hamiltonian of the form 
shown in Eq. (A4). 

The simple procedure for finding the solution is sum¬ 
marized here: 


I. Start with H = 0. Find ujd by requiring |/0 )q and 
\gl)D to be degenerate. 
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2. Find the eigenstates w^)) of the Hamil¬ 

tonian H{n,ujd). For the first step when = 
0 and the eigenstates are degenerate, choose 
|$i,2(f^, Wd)) = (|/0 )d ± \gl)D)/^- 


3. Use Eq. (A6) to calculate w^(r2) := dwd(U)/dU 

4. Set (jJd —y uJd ~\~ and U —^ U -t- AH. 

5. Go to step 2. 

Once the solution is known, the effective coupling g(U) 
is calculated from the eigenenergies Ui^ 2 (U, ujd ) of the two 
eigenstates |$i^ 2 (U,Wd)). Inspection of Eq. (A4) shows 


that these eigenenergies are equal to EoSset ± 9 and there¬ 
fore 


g{n) = 


Ei{n,LUd) - U2(U,Wd) 


Appendix B: Calculation of a± and /3± 

In this appendix we use perturbation theory to cal¬ 
culate the overlaps a±, j3± of the driven Hamiltonian 
eigenstates |$j-) corresponding to the undriven eigen¬ 
states = {\fO)D + \9 ^)d)/V^- The result is then 


used to derive Eq. (24). We retain only terms up to 


second order in the drive strength ft and zeroth order in 
the Jaynes-Cummings coupling g. In this approximation, 
we can treat states with different numbers of photons as 
decoupled and replace the dressed states \ij)D by the cor¬ 
responding bare states \ij). The overlaps a± = ($±|/0) 
and /3± = (4>j- jgl) which we wish to calculate are given 


by 


«± = ^ (i + (/o^'Vo)) 




/3±=±^(l + (5l^'^l5l) 






where \ij ) are the fc-th order corrections to the un¬ 
driven eigenstates \ij). Specifically, 

I/O™) ^ 

-^(1) 

' ' A + a' 

which results in the following expressions for a± and /3±: 


1 


a± = 


/3± = 


1 / y2H/2\ 1 / v^H/2y 

A j ~ 2 yA-a) 
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